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1. Introduction 

Higgs production in association with a jet of high transverse energy with a subsequent 
decay into two photons, pp — > iJ+jet — > 77 + jet, is considered a very promising discovery 
channel for a Higgs boson of intermediate mass (100 GeV < m H < 140 GeV) [1,2]. In fact 
if a high-pT jet is present in the final state, the photons are more energetic than in the 
inclusive channel and the reconstruction of the jet allows for a more precise determination 
of the interaction vertex. Moreover, the presence of the jet offers the advantage of being 
more flexible with respect to choosing suitable acceptance cuts to curb the background. 
These advantages compensate for the loss in the production rate. The analysis presented 
in Refs. [1,2] was done using leading-order perturbative predictions for both the signal [3] 
and for the background [4], although anticipating large radiative corrections which were 
taken into account by using a constant ET NLO = 1.6 factor for both the signal and the 
background processes. 

Since the first analysis was proposed, the radiative corrections to the signal have been 
computed [5]. However, the irreducible pp — * 77 + jet background is still known only at 
leading order, where it proceeds via the sub-processes qq — > 775 and qg — > 77g and is 
dominated by the latter that benefits from the large gluon luminosity. The quark-loop 
mediated gg — > 577 sub-process, which is 0(a|) and thus formally belongs to the next- 
to-next-to-leading order (NNLO) corrections, might have been significant due to the large 
gluon luminosity. It has been computed and found to yield a modest contribution [6, 7] . 
Thus it is not unreasonable to expect that a calculation of the unknown next-to-leading 
order (NLO) corrections provide a reliable quantitative estimate of the irreducible pp — > 
77 + J e t background. In this paper we present the NLO corrections to pp — > 77 + jet. 
Although this background will eventually be measured at the Large Hadron Collider, it is 
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still important to have a better understood theoretical prediction in order to optimize the 
various detector-dependent selection and isolation parameters used for Higgs searches. 

The necessary ingredients for the computations presented here have been known in the 
literature for some time. At NLO accuracy we must consider the tree-level subprocesses 
with two partons and two photons in the final state, namely qq — » 5577 and qq — > QQ77 
together with all the crossing-related subprocesses. The corresponding contributions are 
termed real corrections. We express the matrix elements in terms of amplitudes of fixed 
helicities, compute the square and the sum over helicities numerically. For maximally- 
helicity- violating configurations, the amplitudes were first derived in Ref. [8]. The remain- 
ing helicity configurations were computed in Ref. [9]. In our calculations we used the 
amplitudes presented in Ref. [10], where all helicity amplitudes for processes involving n 
photons and 6 — n partons (n = 1, 2, 3, 4) were re-calculated using the notation of Ref. [11]. 

At NLO we also need the interference between one-loop and tree contributions to the 
Born-level subprocesses qq — > 577 and qg — > 773, which constitute the virtual corrections. 
The one-loop helicity amplitudes are decomposed in terms of primitive amplitudes for the 
two-quark three-gluon one-loop amplitudes [12,13]. The explicit form of the decomposition 
was presented in Ref. [10]. 

In a NLO computation both the real and the virtual corrections are separately diver- 
gent in four dimensions, but their sum is finite for infrared-safe observables. There are sev- 
eral methods to compute this finite correction. We used the dipole subtraction method [14] 
slightly modified for better numerical control [15] as implemented in the NL0JET++ pack- 
age [16]. 

In Sect. I the cross sections for pp — > 77 + jet production at leading order and at 
NLO are described; in Sect. |3] several distributions of eventual phenomenological interest, 
as well as the behaviour of the cross section under renormalization and factorization scale 
variations, are analysed, with and without radiative corrections; in Sect. || we present our 
conclusions. 

2. Cross section formulae 

In this section we spell out the explicit cross section formulas that are used in our computa- 
tion. First we give the leading order production rate, then the rates used in the calculation 
of the radiative corrections: the real and virtual contributions and and certain Born-level 
matrix element expressions needed for the construction of the dipole subtraction terms. 

2.1 The leading order production rate 

Here we summarize what is necessary to compute pp — > 77 + jet at leading order. All the 
relevant parton subprocesses can be obtained, directly or through crossing symmetry, from 
the tree-level scattering amplitude for qq — > 775. The colour decomposition of amplitudes 
involving two quarks, two photons and a gluon is [11] 

4 ree (l 5 ,2 g ;3 9 ;4 7 ,5 7 ) = 2g s e 2 Q 2 q (T^)% Af™^ 2 q ; 3 g ; 4 7 , 5 7 ) , (2.1) 
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where eQ q is the electromagnetic charge of the quark q, g s is the strong coupling and 
T a are the generators of the SU(N C ) algebra in the fundamental representation*. On the 
right-hand side of Eq. fl2.1|) A^ ee is the colour-stripped sub-amplitude which we give for 
fixed helicities of the external particles, with all particles taken as outgoing. Configurations 
with all of the external particles of equal helicity, or all but one, do not contribute at tree 
level. For maximally-helicity-violating configurations, (— , —,+,..., +), the sub-amplitudes 
are [8,11,17] 



^fl+2"-3-4 5 )-i {U){2t)3 {21) (21) 
5 1 9 "' 7,b7j " (12) (2 3) (31) (2 4) (41) (2 5) (51) 



where we make explicit only the helicity of the quark line, and where i is the gluon or photon 
of negative helicity. Eq. ( |2.2| ) defines the only independent sub-amplitude, all of the other 
sub-amplitudes are obtained by relabeling and by use of the discrete symmetries of parity 
inversion and charge conjugation. Parity inversion flips the helicities of all particles. It 
is accomplished by the "complex conjugation" operation (indicated with a superscript >) 
defined such that (ij) <-> [ji], but explicit factors of i are not conjugated to — i. A factor 
of —1 must be included for each pair of quarks participating in the amplitude. Charge 
conjugation changes quarks into anti-quarks without inverting helicities. In addition, there 
is a reflection symmetry in the colour ordering such that 

^-5 ree (l?! 2g! 3 9 ; 4 7 , 5 7 ) = —A^ ee (2 q , l q ; 3 9 ; 4 7 , 5 7 ) . (2-3) 

Using Eqs. fl2.1| ) and ([E2]), the squared amplitude for qq — > 775, summed over colours and 
helicities, is 



V |4 ree (Lj,2 9 ;3 g ;4 7 ,5 7 )| 2 = 8e 4 Q q g 2 s N c C F 



S12 S12 



col, hel 



S12S23S13 S14S24 S15S25 
Sl3S23(sl 3 + S2 3 ) + S14S 2 a(Su + sl 4 ) + Sl 5 S25(Sl 5 + 825) ■ (2.4) 

In order to obtain the production rate for qq — > 775, we cross to the physical region and 
take 1 as the incoming quark and 2 as the incoming antiquark, 

da(qq) = 1 dV 3 \ ^ |^ ree (1,-, 2 q ; 3 g ; 4 7 , 5 7 )| 2 , (2.5) 

c col, hel 

with s = (pi + P2) 2 , and where we make explicit the average over initial colours and 
helicities, and the symmetry factor for the final-state photons. dT^ is the 3-particle phase 
space for the final state, 

m = n jSfiWi (27r)4 6 ' {pi + P2 - P3 - PA - P5) ' (2 - 6) 



*We normalise the fundamental representation matrices by taking tr(T a T 6 ) = Tn8 ab , with Tr = 1. For 
the d quadratic Casimir operators we have CV = Tr(N^ — i)/N c = (Nc — 1)/N C in the fundamental and 
Ca = 2 TrN c = 2 N c in the adjoint representation. 



-3- 



that we generate uniformly. 

The production rates for the subprocesses qg — * (777 and qg — > q-yy are obtained from 
Eq. (2.5) by a suitable relabelling (1 <-> 3 and 2 <-> 3, respectively) of the entries in the 
squared amplitude (|2.4| ), and by replacing the colour average iV 2 with 7V c (iV 2 — 1). 

2.2 The NLO production rate 

At NLO we must consider the tree-level subprocesses with two partons and two photons 
in the final state, namely qq — > ggjj and qq — > and all of the crossing-related 

subprocesses, and the interference between one-loop and tree contributions to qq — > 577, 
and crossing-related subprocesses. Let us first consider the subprocesses with four final- 
state particles. 

For qq — > 5577, the colour decomposition is 

^tree^_^ 2^ 3 S , 4 g ; 5 7 , 6 7 ) = 2<7ge 2 Qg 

(r a 3T a 4) | 2(? ; 3 9 , 4 S ; 5 7 , 6 7 ) + (T"T*)\ Af^lg, 2 q] 4 g , 3 S ; 5 7 , 6 7 )] .(2.7) 

The sub-amplitudes for the maximally- helicity- violating configurations (—,—,+,+,+,+) 
have similar forms as in Eq. ( |2.2[ ) [8,11,17]. Those for the configurations (—,—,—,+,+,+) 
are more complicated, and have been computed in Refs. [9,10]. For each helicity configu- 
ration, the squared amplitude for qq — > 5577, summed over colours, is 

|^ t 6 ree (l 9 -,2 9 ;3 g ,4 s ;5 7 ,6 7 )| 2 = 4e 4 Q^| C F 

col 

x [n c C f [|4 ree (l 5 , 2 q ; 3 g ,4 g ; 5 7 , 6 7 )| 2 + [4 ree (l 5 , 2 9 ; 4 fl , 3 9 ; 5 7 , 6 7 )| 2 ] 

-2Re [4 ree (l 9 -,2 9 ;3 s ,4 g ;5 7 ,6 7 )M t 6 ree (l 5 ,2 (? ;4 s ,3 5 ;5 7 ,6 7 )] } . (2.8) 



As in Sect. 2.1, in order to obtain the production rate for qq — > ^yjgg we cross to the 
physical region and take 1 as the incoming quark and 2 as the incoming antiquark, 

1 AT? 1 1 



da ^ q q^ ligg) = l_ m 1 _L |^ t 6 ree (l 5 ,2 g ;3 5 ,4 9 ;5 7 ,6 7 )| 2 , (2. 

col, hel 



where we made explicit the symmetry factor for the final-state photons and gluons. dV^ is 
the 4-particle phase space for the final state, 

6 ,3 

dV 4 = J[ {27t) s2 Ei ( 2vr ) 4 +P2-P3-P4-P5- Pe) , (2.10) 

which we generate from the three-parton phase space Eq. ( |2,6| ) with a subsequent dipole 
splitting [14]. The production rates for the subprocesses qg — * qg"f~f and qg — > qgj'y and 
99 ~^ 9977 are obtained from Eq. (^1]) by a suitable relabelling of the entries in the squared 
amplitude (|2.8| ) , by changing the symmetry factor to that of photons only, and by a suitable 
replacement of the colour average. 
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The colour decomposition for the subprocesses with two pairs of quarks of different 
flavour, qq — > QQjj, is 

^r(l5>2 g ;3Q,4 Q ;5 7 ,6 7 ) = 2e 2 g 2 s (T<>)| (T«)g < ee 2 q ; 3 Q , 4 Q ; 5 7 , 6 7 ) . (2.11) 



Note that in Eq. (2.11) the electromagnetic charges do not factor out, and thus are 
kept within the sub-amplitudes. The sub-amplitudes for the maximally-helicity-violating 
configurations (— ,—,+,+,+, +) are listed in Ref. [10,11]. Those for the configurations 
(— , — , — , +, +, +) can be found in Ref. [9, 10]. For each helicity configuration, the squared 
amplitude for qq — ► QQ77, summed over colours, is 

l4 ree (l ? - ; 2 9 5 3q, 4q; 5 7 , 6 7 )| 2 = 4e 4 g A s N c C F \A% ee (l q , 2 q ; 3 Q , 4 Q ; 5 7 , 6 7 )| 2 . (2.12) 

col 

The production rate for qq — > QQ77 has the same form as Eq. (|2.9| ), up to the symmetry 
factor which is for photons only, 



da(qq^QQ 11 ) = ^dV i l -^ i £ |^ ree (l,-, 2,; 3q, 4 q ; 5 7 , 6 7 

c col, hel 



(2.13) 



The production rates for qQ — > (7Q77, and for the substitution of one or of both of the 
incoming quarks with antiquarks, are obtained by suitable relabelings of the entries in 
Eqs. ( fn2D and ( gig ). 

For two quark pairs of equal flavour, we must antisymmetrise Eq. Q2.ll ) by subtracting 
the same expression with the colour and momentum labels of the quarks exchanged, 



^6 ree (lg, 2 9 ;3g,4 9 ;5 7 ,6 7 ) = 2e 2 gj 



(T°)24 ree (l 9 -,2 9 ;3 f ,4 g ;5 7 ,6 7 ) 



-(T a )|(T a )K ee (l f ,4 9 ;3 5 ,2 9 ;5 7 ,6 7 ) 



• (2-14) 



Then for each helicity configuration the squared amplitude for qq — » qq'j'j, summed over 
colours, is 



2 l4 ree (l ?J 2,5 3 5 , 4 9 ; 5 7 , 6 7 )| 2 = 4e 4 ^(JV? - 1) 



col 



x l l^6 rce (l<?' 2(?j 3g, 4 g ; 5 7 , 6 7 )| + |-A 6 ree (l g , 4 g ; 3 q , 2 q ; 5 7 , 6 



+^ 2 fc4^Re[4 ree (l g -,2 g ;3 g -,4 g ;5 7 ,6 7 )*4 ree (l 5 ,4 g ;3 5 ,2 g ;5 7 ,6 7 )] \ , (2.15) 



where the delta function 5h 2 h A reminds that the interference term is present only when the 
equal-flavour quarks are indistinguishable, i.e. have the same helicity. The production rate 
for qq — > qqjj has the same form as Eq. ( 2.13] ), where we use Eq. ( |2.15 ). The production 
rates for qq — > qq'j'j, and for the substitution of one or of both of the incoming quarks with 
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antiquarks, are obtained by suitable relabelings of the entries in Eqs. ( ]2.13 ) and ( 2. 15| ) . 
However, note that when the incoming fermions are both quarks or both antiquarks, the 
final-state symmetry factor is 1/4 instead of 1/2. 

When computing the virtual corrections we need the one-loop amplitude for the qq — ► 
77g process. Its colour decomposition is similar to that of the tree-level amplitude given 
in Eq. Q, 

4- loop (l g -,2 g ;3 s ;4 7 ,5 7 ) = 2e 2 g 3 s (T a3 )} 2 A 5:1 (l q , 2 q - 3 9 ; 4 7 , 5 7 ) , (2.16) 

where the one- loop colour sub-amplitude A&i is given in terms of sums of permutations of 
the primitive amplitudes for the qq — ► ggg process [10, 13]. Explicitly, 

^5;i(l(?; 2 g ; 3^; 4 7 , 5 7 ) = — Q q A 5 ^ \l q , 3 g , 2 q , a^g, cr^g) 



E 

Q 2 f 



tr 



L >[~] 



,(2.17) 



where tr Q j is a short-hand for the sum of squared charges of different flavour appearing 

in the quark loop. We only need the primitive amplitudes for the helicity configurations 

L [-] 

of type (—,—,+,+,+)'. The unrenormalized primitive amplitudes A*' 2 can be found 

R Ml 

directly in Ref. [12]. The A 5 must be first rewritten in terms of another set of primitive 
amplitudes (see Ref. [10,12])*. 

The amplitude ( 2.17|) is renormalized by adding the MS ultraviolet counterterm, 



"^5 (!<?> ^<?' ^3' ^7> 5 7 , 

with fa = (11JV C - 2N f )/3, and 



-cr 5 i§^ t 5 ree (l 9 -,2 g ;3 9 ;4 7 ,5 7 ) 



cr 



1 r(l + e)r 2 (l-e) 



(2.18) 



(2.19) 



(4vr) 2 - e T(l - 2c) 

The interference term between the one-loop and the tree amplitudes, summed over colours, 
is 



2 E Re 



A 



l-loop / 



col 



l " (fgi 2^; 3g,; 4 7 , 5 7 )^4g (1^, 2 9 ; 3^; 4 7 , 5 7 ) 

= 8e 4 Q 2 q g A s N 2 C F Re [A 5:1 (l q , 2 q ; 3 g ; 4 7 , 5 7 )4 ree (l 9 -, 2 q ; 3 g ; 4 7 , 5 7 )*] . 
Finally, the one-loop production rate for qq — > (777 is given by 

^l-loop/i n .0 .4 c \Atree(i_ o .a .4 c 

"^•5 V -*-<?> °g> ^7> tJ 7y ,/1 5 I- 1 ?) ■ 6 <7) °gi ^7> u " 



(2.20) 



c coLhel 



2.S 



(2-21) 

^The amplitudes for the configurations with all equal helicities, or all but one, as well as the amplitudes 

for gg — > (777 contribute at one loop, however they do not appear in the NLO production rate since their 

tree level counterparts vanish. 

"•"In Ref. [12] the amplitudes Asa are provided in the dimensional reduction scheme [18-20]. If it is 

wished to have them in conventional dimensional regularization, one must add to Eq. (2.17) the term 

A% = -cr(Nl - l)/(2iY c 2 )A t 5 rco , with c r in Eq. (Eyl). 
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The one-loop production rates for qg — > qj'y and qg — > qjj are obtained by relabeling in 



Eq. (2.20) and by using the suitable colour average. 

In addition to the NLO real and loop matrix elements, for the complete NLO calcula- 
tion we also need (i) a set of independent colour projections of the matrix element squared 
at the Born level, summed over parton polarizations, and (ii) an additional projection of the 
Born level matrix element over the helicity of the external gluon in four dimensions. These 
are required for the construction of the subtraction terms [14]. The Born level process 
involves three coloured partons. In this instance, the colour structure exactly factorises, 
i.e. the colour-correlated squared matrix elements I^W^'^ | 2 are proportional to the Born 
squared matrix element, 

|^ 5 ' fc (l g -,2 <? ;3 s ;4 7 ,5 7 )| 2 = c ik £ |4 ree (l 9 -, 2,; 3 9 ; 4 7 , 5 7 )| 2 , (2.22) 

col, hel 

where Cj& = c^i and c\2 = 1/N C , C13 = C23 = — iV c . The projection of the Born level matrix 
element over the helicity of the external gluon is 

^^(lq, 2 g ; 3+ ; 4 7 , 5 7 )A^ ee (lq, 2 q ; 3 g ;4 7 ,5 7 )* = Se A Q q g 2 s N c CF 

col, hel 

A^ ec {lq , 2 q ; 3jj~; 4 7 , 5 7 )A^ cc {Y^,2 q ; 3 g ; 4 7 , 5 7 )* 
+A i £ ee (l~t, 2 g ; 3+; 4+, 5 7 )A]? ee (l g t, 2 q ; 3 g ; 4+, 5 y )* . (2.23) 

2.3 Checks of the computations 

In order to find the NLO correction to the Born cross sections of infrared-safe observables 
we used the dipole subtraction method [14] with cut phase space for the subtraction term as 
described in [15]. The details of setting up the partonic Monte Carlo programs based upon 
the cross section formulae described in the previous subsections are well documented in 
Refs. [14] and [15]. Here we record the checks we performed to ensure the correctness of our 
prediction for the NLO correction: (i) we compared the Born cross section in Eq. ( |2.5| ) to 
the Born-level prediction of Ref. [6] and found perfect agreement; (ii) we made three inde- 
pendent coding of the relevant NLO squared matrix elements; (iii) we checked numerically 
that the subtraction terms match the singular behaviour of the real correction pointwise 
in randomly chosen Monte Carlo points; (iv) we found that the total NLO correction is 
independent, as it should, of the parameter that controls the cut in the phase space of the 
subtraction term (for details see Ref. [15]); (v) to utilize the process independent feature of 
the dipole method, we used the same program library NL0JET++ that was used for other, 
independently tested NLO computations, only the matrix elements were changed. 

3. Results 

We now turn to the presentation of our numerical predictions. We show cross section values 
at leading order and at NLO for a Large Hadron Collider (LHC) running at 14TeV. The 
values shown at leading order were obtained using the leading order parton distribution 
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functions (p.d.f.'s) and those at NLO accuracy were obtained using the NLO p.d.f.'s of 
the CTEQ6 package [22] (tables cteq611 and cteq6m, respectively). We used the two- loop 
running of the strong coupling at NLO with Aqcd = 226 MeV and one-loop running with 
Aqcd = 165 MeV at leading order. The renormalization and factorization scales were set 
to = = x^fj-o, where for the reference value fiQ we used 



with m 77 the invariant mass of the photon pair. Eq. ([Tl]) reduces to the usual scale choice 
for inclusive photon pair production if Pj c tT vanishes. Our prediction for the 77 + jet 
production cross section is intended for use in the detection of a Higgs boson lighter than 
the top quark, therefore, we assume 5 massless flavours and restrict all cross sections to the 
range of 80 GeV < m 77 < 160 GeV. The electromagnetic coupling is taken at the Thomson 
limit, a em = 1/137. 

Our main goal in this paper is to show the general features of the radiative corrections 
to the cross sections. Therefore, we use a jet reconstruction algorithm and a set of event 
selection cuts, expected to be typical in Higgs searches. In particular, in order to find the 
jet, we use the midpoint cone algorithm [23] with a cone size of R = \J Arj 2 + A(f> 2 = 1, 
with A77 the rapidity interval and A(j) the azimuthal angled Then, we require that both 
the jet and the photons have > 40 GeV and rapidity within \rj\ < 2.5. These are the 
same selection cuts as used in Ref. [6] for computing the gluon initiated 0(a|) corrections. 
Furthermore, we isolate both photons from the partons in a cone of size i? 7 . 

At NLO the isolated photon cross section is not infrared safe. To define an infrared 
safe cross section, one has to allow for some hadronic activity inside the photon isolation 
cone. In a parton level calculation it means that soft partons up to a predefined maximum 
energy are allowed inside the isolation cone. 

The standard way of defining an isolated prompt photon cross section, that matches the 
usual experimental definition, is to allow for transverse hadronic energy inside the photon 
isolation cone up to i^T.max = EP7T1 with typical values of e between 0.1 and 0.5, and 
where p 7 x is taken either to be the photon transverse momentum on an event-by-event 
basis or to correspond to the minimum value in the p 7 T range. In perturbation theory 
this isolation requires the splitting of the cross section into a direct and a fragmentation 
contribution. The precise definition of the two terms, valid to all orders in perturbation 
theory, is given in Ref. [24], where it was shown that a small isolation cone for the photon 
leads to unphysical results in a fixed order computation. For a small cone radius i? 7 , an 
all-order resummation of a s \n.{l/ R 2 ) terms combined with a careful study of the border 
line between perturbative and non-perturbative regions has to be undertaken. To avoid 
such problems, we use _R 7 = 0.4 as the default radius. 

In this work, we deal with only the direct production of two photons, therefore, we use 
a 'smooth' photon-isolation prescription which depends on the production of direct photons 
only and is independent of the fragmentation contribution [25]. This isolation means that 
the energy of the soft parton inside the isolation cone has to converge to zero smoothly 

§ In our NLO computation the midpoint and seedless cone algorithms yield identical cross sections. 
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MO = 4 ( m 77 +PjetTj , 
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p p -^77+ jet, E cm = 14 TeV 
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Figure 1: The invariant mass distribution of the photon pair at LHC energy. The photons and 
the jet are required to have transverse momentum |px| > 40 GeV and lie in the central rapidity 
region of \rj\ < 2.5. The jet is reconstructed according to the midpoint algorithm. 



if the distance in the rj — (p plane between the photon and parton vanishes. Explicitly, 
the amount of hadronic transverse energy Et (which in our NLO partonic computation is 
equal to the transverse momentum of the possible single parton in the isolation cone) in 
all cones of radius r < i? 7 must be less than 

(1 — COS V \ 
i 5~ ' ( 3 - 2 ) 
1 — cos ii 7 / 

where we chose to use n = 1 and e = 0.5 as default values, and we take pyr to be the 
photon transverse momentum on an event-by-event basis. 

In Fig. [l] we plot the invariant mass distribution of the photon pair. Here we see 
the continuum background on which the Higgs signal is expected to manifest itself as a 
narrow resonance in the intermediate-mass range. The dotted (red) line is the leading order 
prediction and the solid (red) one is the differential cross section at NLO accuracy. The 
striking feature of the plot is the rather large correction. The large corrections are partly 
due to the appearance of new subprocesses at NLO as can also be read off the figure. The 
gluon-gluon scattering subprocess begins to appear only at NLO accuracy, and therefore 
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it is effectively leading order. It is shown with a long dashed-dotted (magenta) line: it 
yields a very small contribution. The bulk of the cross section comes from quark-gluon 
scattering both at leading order and at NLO, shown with sparsely-dotted (blue) and long- 
dashed (blue) lines. The quark-quark scattering receives rather large corrections because 
the leading order subprocess can only be a quark-antiquark annihilation process, shown 
with short-dashed (green) line, while at NLO, shown with short dashed-dotted (green) 
line, there are (anti-)quark-(anti-)quark scattering subprocesses. Thus at NLO the parton 
luminosity is sizeably larger. In addition, more dynamic processes are allowed, which 
include t-channel gluon exchange. These contribute to enlarge the cross section in phase 
space regions which are disfavoured at leading order. 

A part of the large radiative corrections is accounted for by the new subprocesses; 
another part is due simply to the enlarged phase space, as can be seen from Fig. ^, where 
the differential distributions in the distance i?j 7 = y\rjj — V-y\ 2 + \4>j ~ <AyF between the 
jet and the photons in the r\-<f> plane are shown, with a selection cut at Rj^ > 0.4. We have 
also studied the distribution of the rapidity separation between the harder photon and the 
jet and we found that it is on average small. Therefore, the peak in the Rj-y h distribution 
near tt suggests that the distribution peaks where the jet and the harder photon are mostly 
back-to-back and central in rapidity. The harder photon is separated from the jet by at 
least Rj lh ~ 2 at leading order, while at NLO they can be as close as permitted by the 
selection cut. On the contrary, in going from leading order to NLO, the distribution for 
the photon with lower transverse momentum, the softer photon, changes in normalization 
but not in shape. 

In the following plots, we require i?j 7 > 1.5, which is also a standard selection cut 
in Higgs searches in the Higgs + jet channel [26, 27]. In fact, on one hand a cut on 
Rj~f 3 > 1.5 affects the leading order and the NLO evaluation in the same way because they 
have the same shape. Thus in that case the cut does not reduce the correction to the m 11 
distribution. On the other hand, a cut on Rjj h > 1.5 cuts the NLO correction, but does 
not cut the leading order, so the correction is reduced. Nevertheless, the reduction is less 
then 10%, thus cutting on Rj^ at Rj~ > 1.5 reduces the NLO correction to the invariant 
mass distribution of the photon pair, but only marginally. 

In Fig. ||, we plot the differential distribution of the transverse momentum of the 
photon pair, p 77 T = Ip 7i ,t +P7 2 ,tIi with a cut at p^l^ = 40 GeV. At leading order the 
jet recoils against the photon pair and the respective jet and photon pair distributions 
are identical. At NLO the extra parton radiation opens the part of the phase space with 
p 77 T < Pj^tT- ^he double peak around 40 GeV is an artifact of the fixed-order computation, 
similar to the NLO prediction at C = 0.75 for the C-parameter distribution in electron- 
positron annihilation. The fixed-order calculation is known to be unreliable in the vicinity 
of the threshold, where an all-order resummation is necessary [28], which would result 
in a structure, called Sudakov shoulder, that is continuous and smooth at p 77 T = Pj^tT - 
Without the resummation, we must introduce a cut, p 77 ,T > 40 GeV to avoid regions in 
the phase space where the fixed-order prediction is not reliable. In the following plots, in 
addition to the default cuts, we also require p 77 ,T > 40 GeV. 

The magnitude of the NLO corrections is also heavily influenced by the photon-isolation 



- 10 - 



2500 



2000 — 



^ 1500 I— 



b 



1000 — 



500 — 







p p ^77+ jet, E cm = 14 TeV 

"i — rp — i — | — i — rp — i — | — i — rp — i — | — i — rp — i — | — i — rp — r 

NLO, harder photon 



LO, harder photon 
NLO, softer photon 
LO, softer photon 



J 



harder photon 



softer photon ,_r 




eJJj I I I Jj i ■ i Jj I I I ill - T: 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



"jet 7 



Figure 2: The distributions in the distances between the jet and the photons in the rj-cj) plane. 
The same selection cuts are used as in Fig. 1. 



cut parameters. In Figs. |4|-|5] we show the effect of scanning the parameters i? 7 and e in 
Eq. ( |3.2| ) over the values of i? 7 = 0.4, 1 and e = 0.1, 0.5, 1. We see that the leading order 
predictions are independent of these parameters^, but the NLO ones depend strongly on 
the photon isolation. Firstly, we note that the smaller i? 7 the larger the NLO corrections. 
In addition, the larger e the larger the NLO correction, with the effect being larger if i? 7 is 
larger. This behaviour is in agreement with Eq. (3^2), according to which smaller and 
larger e imply larger amounts of soft-parton energy that is allowed inside the cone. 

Another remarkable feature of Figs. ^45| is that with the applied cuts, the two-photon 
plus jet background for the search of a Higgs boson with mass in the 120-140 GeV range is 
rather flat, therefore, well measurable from the sidebands around the hypothetical Higgs 
signal. This feature is very different from the shape of the background to the inclusive 
pp — > H — > 77 channel, which is steeply falling. 

In order to assess the stability of the predictions against scale variations, we show 
the cross section in a 3 GeV bin around m 77 = 120 GeV, that is the background for a 
hypothetical Higgs signal for a Higgs particle of mass 120 GeV. Fig. | shows the cross 

^The leading-order cross section does not depend on the chosen values of _R 7 because we have required 
Rj~t > 1-5 and the jet momentum is the same as that of the only parton in the final state. 
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Figure 3: Transverse momentum distribution of the photon pair. The same selection cuts are used 
as in Fig. 1. 



section for two sets of photon-isolation parameters. We show the scale variations for varying 
the renormalization and factorization scales separately, keeping the other scale fixed, as 
well as varying them simultaneously. The lower three curves represent the leading order 
predictions. At leading order the dependence on the renormalization and factorization 
scales is rather small, especially when the two scales are set equal (densely dotted line). 
Observing the predictions we conclude that the scale dependence at leading order does 
not represent the uncertainty of the predictions due to the unknown higher orders. The 
inclusion of the radiative corrections results mainly in the substantial increase of the cross 
section. For our default isolation {R~ = 0.4 and e = 0.5) the NLO corrections at = 1 
are about 120 % of the leading order prediction. In addition, the scale dependence is not 
reduced by the inclusion of the radiative corrections. If we require more stringent photon 
isolation cuts, then we find smaller corrections and a more stable prediction. For instance, 
in Fig. ^ we also show the scale dependence of the cross section obtained with i? 7 = 1 and 
e = 0.1. We find that the cross section at NLO is about 40 % larger than the leading order 
prediction, and in this case the scale dependence at NLO is reduced as compared to the 
one at leading order accuracy. 
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Figure 4: Dependence of the invariant mass 
distribution of the photon pair on the photon 
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tion to the selection cuts used in Fig. 1, we 
also require Rj 7 > 1.5 and p 77 ,t > 40GeV. 



pp^77+jet,E cra = 14TeV 




110 120 130 
M,,. [GeV] 



Figure 5: Dependence of the invariant mass 
distribution of the photon pair on the photon 
isolation parameter e for i? 7 = 1. The same 
selection cuts are used as in Fig. 4. 



4. Conclusions 

In this paper we computed the QCD radiative corrections to the pp — > 77 + jet process 
that constitutes part of the irreducible background to the pp — ► H + jet — » 77 + jet 
discovery channel of an intermediate- mass Higgs boson at the LHC. We used a smooth 
photon isolation that is infrared safe to all orders in perturbation theory and independent 
of the photon fragmentation into hadrons. The predictions were made with a partonic 
Monte Carlo program that can be utilized for a detailed study of the signal significance 
for the intermediate-mass Higgs boson discovery in the Higgs + jet production channel 
at the LHC. We found large radiative corrections, however they are rather sensitive to 
the selection cuts and photon isolation parameters. Choosing a mild photon isolation, i.e. 
a small isolation cone radius = 0.4 with relatively large hadronic activity allowed in 
the cone results in more than 100 % correction with as large residual scale dependence 
at NLO as at leading order. Making the photon isolation more stringent, for instance 
increasing the cone radius to i? 7 = 1 and decreasing the hadronic activity in the cone 
reduces both the magnitude of the radiative corrections as well as the dependence on the 
renormalization and factorization scales. This result shows that a constant X NLO = 1.6 
factor, as mentioned in the Introduction, is certainly not appropriate for taking into account 
the radiative corrections to the irreducible background of the pp — » H + jet — ► 77 + jet 
discovery channel at the LHC. 
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